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Several lifetimes in 23 Mg have been determined for the first time using the Doppler-shift atten¬ 
uation method. A Monte Carlo simulation code has been written to model the 7 -ray line shape. 

An upper limit of r < 12 fs at the 95% C.L. has been obtained for the astrophysically important 
7787 keV state. 


I. INTRODUCTION 

The radionuclide 22 Na is thought to be produced in 
significant amounts in ONe-type classical novae [l|. It 
decays by p + v e emission to the first excited state in 22 Ne 
with a mean lifetime of 3.75 years. The daughter nucleus 
quickly relaxes to the ground state by emitting a 7 ray 
with a characteristic energy of 1275 keV. The lifetime of 
3.75 years is sufficiently short to ensure that the ejected 
22 Na is still close to the nova by the time the decay oc¬ 
curs, making it possible to correlate the 7-ray intensity 
with other observational properties, yet sufficiently long 
for the 22 Na to survive beyond the opaque phase of the 
explosion, allowing the 7-ray to be observed. These for¬ 
tunate circumstances make 22 Na the most promising can¬ 
didate for detection of nuclear 7-ray lines from novae Q ■ 
Nova models predict 22 Na 7-ray fluxes 1 order of mag¬ 
nitude below the observational upper limit Q. While 
detection does not appear to be imminent, it seems rea¬ 
sonable to assume that it will happen within the fore¬ 
seeable future [H, and consequently reliable estimates of 
22 Na yields in novae are desirable. Nova production of 
22 Na may also be relevant for explaining the nonstandard 
abundance of 22 Ne observed in Ne-E meteorites [r>|. 

Under the conditions found in novae 22 Na is mainly 
destroyed by proton capture. Two direct measurements 
of the proton-capture cross section have been reported. 
Such measurements are very challenging because the tar- 
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get material is radioactive. The first measurement was 
performed at Ruhr-Universitat Bochum, Germany, in two 
attempts with results published in 1990 @| and 1996 @- 
The second measurement was performed at the Univer¬ 
sity of Washington in Seattle with results published in 
2010 [§j] and 2011 Both studies conclude that at 
peak nova temperatures the reaction is dominated by a 
resonance at a proton energy of E p = 213 keV, corre¬ 
sponding to an excitation energy of E x = 7787 keV in 
the compound nucleus 23 Mg. (At the very highest nova 
temperatures a resonance at 288 keV also makes a sig¬ 
nificant contribution.) However, the resonance strengths 
reported for the 213 keV resonance differ substantially: 


wy = 1.8 ± 0.7 meV (Bochum) 
cj 7 = 5.7+j/g meV (Seattle) 


This leads to an uncertainty of about a factor of 2 in the 
amount of 22 Na ejected from novae [10. 

Since p+ 22 Na and 7+ 23 Mg are the only open channels, 
the resonance strength is given by 


LO 7 = 


2J+ 1 


(2 Jp + 1)(2 JNa + 1) T 


-B p ( 1 - B p ) 


(1) 


where J p = 1/2 and JNa = 3 are the spins of the proton 
and 22 Na, respectively, J is the spin of the resonance, r is 
the mean lifetime, and B p is the branching ratio for pro¬ 
ton decay. Thus a measurement of J, t, and B p provides 
an indirect determination of the resonance strength. The 
spin-parity and mean lifetime of the 7787 keV state have 
been determined to be J n = 7/2^+) and r = 10 ±3 fs, re¬ 
spectively, in a measurement of the fusion-evaporation re¬ 
action i2 c( i2 c,n) 23 M g mm. The branching ratio for 
proton decay has been determined to be B „ = 0.037±0.07 
in a measurement of the f3p decay of 23 A1 mi, which also 
provides compelling evidence for a J 7r = (7/2) + assign- 
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FIG. 1. Decay scheme of the astrophysically important 
7787 keV state. The levels are labeled by their excitation 
energy in keV. 


ment. Thus the indirect approach gives 

w 7 = 1.3 ± 0.5 meV (Indirect, J = 7/2) 

with the uncertainty obtained by adding the uncertain¬ 
ties on r (30%) and B p (19%) in quadrature. This value 
agrees with the old and less precise result of the Bochum 
group, but disagrees with the more recent and more pre¬ 
cise result of the Seattle group. To resolve the disagree¬ 
ment we set out to remeasure r. 

We note that recently a different spin-parity assign¬ 
ment for the 7787 keV state has been proposed by Tri- 
pathi et al. 0 Based on a comparison of the exper¬ 
imentally determined /3-decay strength of the 7787 keV 
state to the /3-decay strength predicted by a shell-model 
calculation, Tripathi et al. argue that J 71 ' = 5/2 + is a 
more plausible assignment, which gives 

wy = 1.0 ± 0.4 meV (Indirect, J = 5/2) 

aggravating the discrepancy with respect to the Seattle 
measurement. 

In a recent measurement of the /3y spectrum of 
23 Al [l5|, iH| the 7787 keV state is observed to decay to 
the 451 keV and 2052 keV states with relative intensities 
of 81(4)% and 19(4)% and 7 -ray energies of 7335.2(6) keV 
and 5735.4(7) keV. The excitation energy of the 7787 keV 
state is determined to be E x = 7787.2(6) keV, which 
is the most precise value to date and the value that 
we adopt in the present study. While originally only 
published in a PhD thesis pjj, this value has since 
been quoted in a peer-reviewed publication fl3l |. We 
further note that the value agrees with the value of 
E x = 7785.7 ± 1.1 keV obtained in the 12 C( 12 C, n) 23 Mg 
experiment fill. fl2|. where only the 7787 keV —> 451 keV 
transition was observed. (The value given in Ref. 0 is 
E x = 7784.6 ± 1.1 keV, but this value appears not to have 
been corrected for the nuclear recoil.) A simplified level 
scheme is shown in Fig. [T] 

II. EXPERIMENTAL TECHNIQUE 

We employ the Doppler-shift attenuation method 
(DSAM) []j|. The excited nucleus is produced at 


high speed and slows down in a dense medium. If the 
lifetime is sufficiently short, the nucleus will decay before 
it comes to rest, and hence the energy of the emitted 
y-ray will be Doppler-shifted according to, 


E 7 — £ 7 ,o 


(1 — /3 2 ) 1/2 

1 — /3 cos # 7 


( 2 ) 


where £ 7,0 is the 7 -ray energy at rest, /3 is the instan¬ 
taneous speed of the nucleus in units of c, and 0 7 is the 
angle of emission of the 7 ray relative to the trajectory 
of the nucleus. A shorter lifetime leaves less time for the 
nucleus to slow down resulting in a larger Doppler-shift 
and vice versa. Consequently, the 7 -ray line shape is sen¬ 
sitive to the lifetime of the excited nucleus. In addition 
to the lifetime, several experimental effects also influence 
the 7 -ray line shape. These are discussed in Section IIVI 
For very short-lived excited states the average speed re¬ 
duction is roughly 5/3 ss 0.32 x 10 -3 x (r/A) x (d£/da;), 
where r is the mean lifetime in units of fs, A is the 
mass number of the excited recoil, and d£/da; is the 
stopping power in units of MeV/pm. Three important 
observations can be made from this formula: First, a 
dense medium with a high stopping power gives the 
highest sensitivity. This makes Au an ideal stopping 
medium. Second, the stopping power must be known 
to extract the lifetime. At the speeds relevant to the 
present study ( [3 = 0.059-0.082) the energy loss of Mg 
ions in Au is dominated by electronic stopping, and the 
stopping powers have been determined experimentally 
by Forster et al. 0 to a precision of ±4%. Approxi¬ 
mately, dE/dx ~ 8.5 MeV/pm. Third, lifetimes shorter 
than 1 fs will be very difficult to measure even in an 
idealized experiment. Assuming r = 1 fs one obtains 
5/3 « 1.2 x 10 -4 , which implies a shift of only 0.9 keV for 
the 7787 keV —>-451keV transition. 


III. SETUP AND PROCEDURE 

The experiment was carried out at the ISAC-II facility 
of TRIUMF. A schematic drawing of the setup is shown 
in Fig. [2] A similar setup has been successfully employed 
in previous experiments at TRIUMF to measure lifetimes 
of excited states in 19 Ne [2(| 0| and 15 O (22|. The 
24 Mg beam is accelerated to 3.125 A MeV corresponding 
to (3 = 0.082 and directed onto a 3 He-implanted Au foil. 
Excited states in 23 Mg are populated via 3 He( 24 Mg,a) 
with a Q-value of 4.05 MeV. Kinematic curves are shown 
in Fig. [31 Because the experiment is performed in inverse 
kinematics, the recoils are strongly focused in the forward 
direction with speeds in the range (3 = 0.059-0.065. The 
a particles are detected using a A E-E telescope placed 
downstream of the target, consisting of two Si surface- 
barrier detectors with an active area of 150 mm 2 and 
thicknesses of 100 pm and 500 pm. The target is suffi¬ 
ciently thick to stop the beam and the recoils, thus pre¬ 
venting damage of the detectors, while thin enough to 
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FIG. 2. Schematic drawing of the experimental setup. The setup also includes a cooling system and a collimator not shown 
on the drawing. The collimator consists of two circular apertures with diameters of 2.5 and 3.0 mm placed at distances of 73 
and 49 mm, respectively, upstream from the target. e 7 is the energy- and angle-dependent detection efficiency of the HPGe 
detector, here shown for a fixed energy of I? 7 = 8 MeV. 
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FIG. 3. Target composition and 3 He depth profile predicted by SRIM [23] . The surface contamination resulted from conden¬ 
sation of residual gas and water vapor. 


allow transmission of the a particles. A 4 mm diameter 
circular aperture placed at a distance of 14.2 mm down¬ 
stream from the target limits the angular acceptance of 
the A E-E telescope to 0-8°. A 7 -ray detector consist¬ 
ing of four separate HPGe crystals is situated 7.83 cm 
downstream of the target and covers angles up to ps 32°. 
The placement of the detectors at 0° serves to maximize 
the Doppler shift. By gating on the a particles we sup¬ 
press competing channels such as n, p, 2 p, d, n + p , 3 He, 
d+p, and n + 2p, resulting in a significantly cleaner 7 -ray 
spectrum. By gating on the a particles we also reduce 
the spread in recoil velocities (magnitude and angle) and 
thus the kinematic broadening of the 7 -ray lines. 

The target was prepared at the Universite de Montreal 
by directing a 30 keV 3 He beam onto 22.6 ± 0.2 /mi thick 
Au foil. In this way the 3 He ions were implanted at shal¬ 
low depth with a number density of 6.6 x 10 17 cm -2 . The 
depth profile predicted by SRIM [23| is shown in Fig. [3] 
It peaks at 0.07 /im, but is quite broad with a tail that 
extends 0.15 /mi into the Au foil. The thickness of the 
Au foil was determined using a precision mass balance. 


The integrated beam time was 5 days. Typical beam 
intensities were 1-2 x 10 1Q ions/s. The target tempera¬ 
ture was kept at —80° using a LN 2 cooling system, see 
Refs. [ 20 M 22 ] for details. This is necessary to prevent loss 
of 3 He due to beam-induced heating of the target. The 
elastic scattering rate was used for on-line monitoring of 
the 3 He content of the target. No indication of 3 He loss 
was seen. 

During the initial pump-down residual gas and water 
vapor was inadvertently condensed on the surface of the 
cooled target due to poor vacuum. The thickness of the 
condensation layer has been estimated to be 0.2 pm of 
water equivalent based on the observed energy of the 
main 3 He and a-particle groups. As a result the beam 
energy was degraded by 0.25 MeV. Reactions between 
the beam and the surface contaminants contribute to the 
background in the 7 -ray spectrum, but serendipitously 
also provide us with two unshifted 7 -ray lines at high 
energy, namely, the 6877.88(8) keV line from 28 Si and 
its single-escape peak, which enable us to determine the 
energy calibration and the intrinsic detector resolution 
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FIG. 4. Kinematic curves for the 8 He( 24 Mg, a) reaction at 
a beam energy of 3.125 A MeV. Hatched regions indicate de¬ 
tected a particles and associated 23 Mg* recoils. 

very precisely in the energy region of interest. (The 7 - 
ray energy given here differs slightly from the 7 -ray en¬ 
ergy of 6877.0 keV given in the NDS evaluation [2dl] . The 
7 -ray energy given here was calculated from the excita¬ 
tion energies given in the NDS evaluation, corrected for 
the nuclear recoil energy.) The 6877.88 keV line and its 
single-escape peak are shown in Fig. [5] The uncertainty 
on the energy calibration is ±1 keV at 7.8 MeV. Gain 
variations during the experiment were below ±0.5 keV. 
The yield of the 6877.88 keV line relative to the elastic 
scattering rate was constant during the experiment, con¬ 
firming that the condensation took place only during the 
initial pump-down of the chamber. 

IV. LINE SHAPE 

Several experimental effects contribute to the 7 -ray 
line shape. In order to extract reliable lifetimes these 
effects must be accurately and precisely modeled. Be¬ 
low, we list the effects and estimate their contribution 
to the relative spread in 7 -ray energy, 8E 1 /E 7 , at E 1 = 
7.8 MeV. The estimates are rather crude and merely serve 



FIG. 5. The unshifted 7 -ray line due to the 6879 keV —► g.s. 
transition in 28 Si. 

to provide an idea of the relative importance of the var¬ 
ious effects. The quantitative treatment is discussed in 
Section [V] 

Intrinsic resolution: The intrinsic resolution of the 7 - 
ray detector is reflected in the line shape of un¬ 
shifted peaks. At the highest 7 -ray energies the 
line shape is found to be well described by a Gaus¬ 
sian; at the lowest 7 -ray energies two additional 
components, a skewed Gaussian and a smoothed 
step function, must be included to obtain an ade¬ 
quate fit, see Ref. [22j for details. Fits to the un¬ 
shifted 7 -ray lines due to the 279 keV -± g.s. tran¬ 
sition in 197 Au (r = 26.8 ps) and the 6879 keV —> 
g.s. transition in 28 Si (r = 2.7 ps) are shown in 
Fig. [Gj The FWHM of the Gaussian component 
is found to obey the relation, FWHM = 3.3 keV x 
{E 1 / MeV ) 0 52 . This gives 8E J /E 1 « 1.2 x 1CT 3 at 
E 1 = 7.8 MeV. 

Angular acceptance: The finite angular acceptance of 
the detectors results in kinematic broadening of 
the Doppler-shifted 7 -ray lines. Accurate model¬ 
ing of this effect requires knowledge of the detec¬ 
tor geometry and the angular dependence of the 
7 -detection efficiency. In principle knowledge of 
the angular distribution of the a particles and the 
a -7 angular-correlation function is also required, 
but their effect on the line shape is rather small 
due to the limited angular acceptance. The fi¬ 
nite acceptance of the A E-E telescope causes a 
spread of 8E 1 /E 1 « 0.6 x 10 -3 , while the finite 
acceptance of the 7 -ray detector causes a spread of 
5E 1 /E 1 « 0.6 x 10" 2 . 

Beam emittance: The finite emittance of the beam 
also contributes to the line shape. The energy 
spread translates directly into a spread in speed and 
thus a spread in the Doppler shift. In the present 
experiment the energy spread of the beam was 
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FIG. 6 . (Color online) Fits to the unshifted 7 -ray lines due to the 279 keV —t g.s. transition in 191 Au and the 6879 keV -7- g.s. 
transition in 28 Si. 


0.2% (FWHM), resulting in a spread of SE^/E-y « 
0.6 x 10~ 4 . The divergence and transverse size of 
the beam influences the line shape in a subtler way 
by allowing the detection of a particles with scat¬ 
tering angles larger than the 8 ° acceptance of the 
circular aperture in front of the A E-E telescope. 
The divergence of the beam can safely be neglected. 
The transverse size is constrained by the 2.5 mm 
diameter circular aperture placed at a distance of 
73 mm upstream from the target. Visual inspection 
of the target at the end of the experiment revealed a 
burn mark with a diameter of 2-3 mm, consistent 
with the size of the aperture. A beam diameter 
of 2.5 mm thus appears realistic, and we conclude 
that a particles with scattering angles as large as 
13° reached the detector. We estimate the spread 
caused by this effect to be dE^/E^ ps 0.9 x 10 -3 . 

3 He depth profile: Since a reaction can take place at 
any depth inside the 3 He-implanted layer, the 
beam-energy resolution is fundamentally limited by 
the energy loss inside this layer. The maximum en¬ 
ergy loss, which occurs if the reaction takes place at 
the bottom of the 3 He-implanted layer near a depth 
of 0.15 /irn, is 1.0 MeV. This results in a spread of 
5E 1 /E 1 « 0.4 x 10" 3 . 

Multiple scattering: Angular deflection due to multi¬ 
ple scattering in the target influences the line shape 
in the same way as the finite beam emittance by al¬ 
lowing the detection of a particles with scattering 
angles larger than the 8 ° acceptance. The half¬ 
width of the angular distribution ( 07/2 in Eq. Q) 
may be estimated from Eq. (20) in Ref. [25], though 
the stated range of applicability (20-90 A MeV) is 
somewhat above the range of speeds probed in the 
present study (5-10 A MeV for the a particles and 


1-2 A MeV for the 23 Mg recoils). For the a par¬ 
ticles the half-width after passage through the Au 
foil is found to vary between 1.8° and 3.6°, the 
precise value depending on the energy, resulting in 
a typical spread of 5E 1 /E 1 m 1.0 x 10 -3 . The 
23 Mg recoils are also deflected, but less so: the 
half-width after a distance of 0.36 /irn in the Au 
foil (the distance covered in approximately 20 fs) 
is found to vary between 0 . 6 ° and 1 . 2 °, the precise 
value again depending on the energy. The lateral 
deflection caused by multiple scattering can safely 
be neglected. 

Energy-loss straggling: The energy-loss straggling 
may be estimated using the semi-empirical formula 
of Ref. [3] • Using the same numbers as above 
(40 MeV, 0.36 pm) the energy loss of the 23 Mg re¬ 
coils is determined to be 3.4 MeV and the straggling 
to be 0.23 MeV (FWHM). This results in a spread 
of 5E 1 /E 1 « 1.8 x 1CT 4 . 

The experimental effects contributing to the 7 -ray line 

shape are summarized in Table [T] 


V. MONTE CARLO SIMULATION 

A quantitative treatment of the experimental effects 
discussed in the preceding section is most easily accom¬ 
plished with a Monte Carlo simulation. A simulation pro¬ 
gram has been written in the language C-|—h The pro¬ 
gram makes use of functionalities provided by the data 
analysis framework ROOT [27| . A flowchart of the pro¬ 
gram is shown in Fig. 0 The program may be obtained 
at https://github.com/oliskir/SimDSAM.git 

We begin by sampling the properties of the beam par¬ 
ticle, he., its energy and transverse displacement from 
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FIG. 7. Flowchart of the Monte Carlo simulation program. See text for details. 


the beam axis, the depth of the implanted 3 He atom in 
the target, the emission angles of the a particle in the 
c.m. frame, the emission angles of the 7 ray in the rest 
frame of the 23 Mg* recoil, and the lifetime of the excited 
state (t). For the beam we assume a Gaussian energy dis¬ 
tribution with a spread of 0.2% (FWHM) and a uniform 
circular transverse profile with a diameter of 2.5 mm. For 
the 3 He implantation depth we assume the profile deter¬ 
mined by SRIM [231] shown in Fig. [3] For the emission 
angles we generally assume isotropic distributions. The 
error introduced by this simplifying assumption has been 
estimated by trying a few different realistic anisotropic 
distributions and has been found to be negligible com¬ 


pared to other errors for all the transitions considered, 
except the 2908 keV —> g.s. transition where it gives an 
error of ±2 fs. Finally, the lifetime is drawn from the 
standard exponential-decay distribution. 

Next, we propagate the beam particle in the target 
until it reaches the specified depth. This involves three 
tasks: First, the energy loss is calculated by numerical 
integration of the stopping powers (cf. Section IVT1) assum¬ 
ing that the particle travels along a straight path. (This 
assumption implies a slight underestimate of the energy 
loss which, however, is negligible at the energies relevant 
to the present study.) Second, the angular deflection due 
to multiple scattering, a, is drawn from the distribution 
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TABLE I. Experimental effects contributing to the 7 -ray 
line shape and their relative importance for the 7786 keV —>• 
451 keV transition in 23 Mg compared to the shift caused by a 
lifetime of 1 fs. 


Experimental effect 

5E 1 /E 1 (% 0 ) 

Intrinsic resolution 

1.3 

A E-E acceptance 

0.6 

HPGe acceptance 

6 

Beam-energy spread 

0.06 

Beam transverse size 

0.9 

3 He depth profile 

0.4 

Multiple scattering 

IIP 

Energy-loss straggling 

0.18^ 

Shift caused by r = 1 fs 

0.12 


a Assuming r = 20 fs. For shorter lifetimes the effect is smaller. 


(Eq. ( 8 ) in Ref. @) 

/ a 2 \ 

/(a) oc sin (a) exp —In 2 —— , (3) 

\ a i/2 / 

where 07/2 is the half-width given by Eq. (20) in Ref. [25:). 
Note that for small deflections sin (a) ss a, implying a 
“universal” distribution of the form /(it) oc wexp(—it 2 ) 
with u = Vln 2 0 / 07 / 2 . Third, the energy-loss straggling 
is sampled from a Gaussian distribution with the stan¬ 
dard deviation given by the semi-empirical formula of 
Ref. & 

We determine the momenta of the a particle and the 
excited 23 Mg* recoil in the centre of mass (c.m.) frame, 
using the already sampled emission angles, and boost to 
the laboratory (lab) frame, always using relativistic kine¬ 
matics. We then propagate the a particle and the 23 Mg* 
recoil in the target, the former until it exits, the latter 
until it decays, i.e., until t = f 'yds/v, where v = c/3 is 
the instantaneous speed of the recoil, 7 = (1 — / 3 2 ) -1 / 2 , 
and ds is an infinitesimal step along the trajectory of the 
recoil. We then determine the momentum of the emitted 
7 ray, again using the already sampled emission angles, 
and boost to the lab frame. 

Finally, we check if the a particle and the 7 ray strike 
the detectors and sample the intrinsic resolution of each 
detector. For the 7 -ray detector we also use the energy- 
and angle-dependent detection efficie ncy e 1 ( E 7 , 9 1 ), de¬ 
termined by a GEANT4 simulation [2|| and shown in 
Fig.[ 2 ]for E 7 = 8 MeV. (In practice, simulations were per¬ 
formed for E 1 = 0.1, 0.5, 1, 2, 3 ..., 10 MeV and the effi¬ 
ciency at intermediate energies was determined by cubic- 
spline interpolation.) 

VI. STOPPING POWERS 

Inside the 3 He-implanted layer the stopping power is 
different from that in pure Au. The presence of the 3 He 
atoms acts to increase the stopping power, but also causes 


the region to swell which lowers the density and thus 
the stopping power. This must be taken into account 
when the energy loss is calculated. Letting c( x) denote 
the local relative concentration of 3 He at depth x in the 
target, the local swelling is given by /S.V/V = Ac{x) 
where A = 0.75 ± 0.25 [29] . We take c(x) to be the 3 He 
depth profile predicted by SRIM [23] normalized to the 
total number of implanted 3 He atoms, 6.6 x 10 17 cm -2 . 
As seen in Fig. [3] the concentration peaks at a depth of 
xo = 0.07 fin 1 , reaching a maximum value of c(x 0 ) = 1.05. 
The stopping power may then be calculated as 

dg = p Au r (dE_\ + 3 /d£\ 1 

dx l + ^4c(x) L \pdxJ Au 197 ; VPdccy 3He 

( 4 ) 

where pAu = 19.30 g/cm 3 is the density of gold and 
(dE/pdx)Au and (dE/pda;)3 Ho are the stopping powers 
in Au and 3 He in units of MeV/(g/cm 2 ). For Au we use 
the experimental stopping powers of Ref. 0 , measured 
to a precision of ±4%. For 3 He we use SRIM stopping 
powers [23). In both cases we interpolate between the 
tabulated values using cubic-spline interpolation. 


VII. DATA ANALYSIS 

As shown in Fig. [5) the differential energy-loss data 
from the A E-E telescope allows us to cleanly separate 
a particles from other light ejectiles, such as p, d , and 
3 He. The reaction channel of interest, ( 3 He, a ), produces 
a particles with a well-defined energy at 0 ° degrees, pro¬ 
viding unambiguous determination of the excitation en¬ 
ergy of the 23 Mg* recoil. Experimental effects smear out 
the measured energy, resulting in an a-particle energy 
resolution of 1 MeV (FWHM), which translates into an 
excitation-energy resolution of about 0.5 MeV (FWHM). 
A standard a source, consisting of 239 Pu, 241 Am, and 
244 Cm, was used to calibrate the Si detectors. 

By gating on the Doppler-shifted 7 -ray line of the 
2908 keV —> g.s. transition in 23 Mg, one obtains the a- 
particle spectrum drawn by the solid line in Fig. [3] (a), 
where a clear peak is seen at E a ss 32 MeV. When the 
energy loss of the a particles in the Au foil is taken into 
account (» 4 MeV), the peak energy is found to be con¬ 
sistent with the excitation energy of the 2908 keV state, 
cf. Fig. [I] The a-particle spectrum obtained by gating on 
the Doppler-shifted 7 -ray line of the 7787 keV - 7 - 451 keV 
transition is shown by the solid line in Fig. [5] (b). Again 
a clear peak is seen, this time at E a ss 21 MeV, an 
energy that is consistent with the excitation energy of 
the 7787 keV state. In contrast, no narrow peak is 
seen in the a-particle spectrum obtained by gating on 
the unshifted 7 -ray line of the 6878 keV —1 g.s. transi¬ 
tion in 28 Si. This strongly suggests that the 28 Si* re¬ 
coil is produced in a fusion-evaporation reaction, such as 
24 Mg+ 12 C —l 2a + 28 Si* or 24 Mg+ 16 0 -a 3a + 28 Si*. The 
7 -gated a-particle spectrum may thus provide important 
clues to the origin of unidentified 7 -ray lines. 
















E (MeV) 


FIG. 8 . (Color online) A E-E plot used to identify the light 
reaction products. Each banana-shaped band represents a 
reaction product or a coincidence of two reaction products. 
From top to bottom: a + a, p + a, a, 3 He, p + d, p + p, d, 
and p. 



FIG. 9. a-Particle spectra obtained by gating on the 
2908 keV — y g.s. and 7787 keV —> 451 keV transitions in 23 Mg 
and the 6878 keV —y g.s. transition in 28 Si. The background 
spectrum shown by the dashed line in (a) has been obtained 
by gating on an interval of 7 -ray energies just above the 7 -ray 
peak. 


By inverting the procedure, we may generate ct-gated 
7 -ray spectra with rather tight cuts on the excitation 
energy in 23 Mg. Thus, feeding from higher-lying states 
is efficiently suppressed, which simplifies the line-shape 
analysis considerably. The a-gated 7 -ray spectra of the 
2908 keV and 7787 keV states are shown in Fig. [T3] and 
Fig. I20I17I18I19I20I The gate on the a-particle energy 
was 2.4 MeV wide in both cases. Simulated line shapes 
have been fitted to the experimental data for three differ¬ 
ent mean lifetimes, using the method of maximum like¬ 
lihood [§3|. The only free parameters are the normal¬ 
ization, two or three parameters to describe the back¬ 
ground, and an energy shift, which is allowed to vary 
within ±(<7t a b + 0 cai) where <Tt a b is the uncertainty on 
the tabulated energy (typically la) and cr ca i is the la 
uncertainty on the energy calibration. The background 
was modeled as a constant plus a smoothed step func¬ 
tion in the case of the 2908 keV —> g.s. transition, and 
as linear function plus a smoothed step function in the 
case of the 7787 keV —y 451 keV transition. The y 2 /dof is 
given for each fit (the number enclosed in parentheses). 
The mean lifetime of the 2908 keV state is determined 
to be 15 ± 3 fs, while an upper limit of 12 fs is obtained 
for the 7787 keV state at the 95% C.L. The overall er¬ 
ror of ±3 fs on the mean lifetime of the 2908 keV state 
includes a systematic uncertainty of ±2 fs from the an¬ 
gular distribution (cf. Section [V]) , which has been added 
in quadrature to the statistical uncertainty from the fit. 
We note that the fit is not too sensitive to the constraints 
imposed on the energy shift. If the constraints are com¬ 
pletely removed, the preferred value for the mean lifetime 
of the 2908 keV state is reduced slightly to 12 fs, and the 
upper limit on the mean lifetime of the 7787 keV state is 
increased slightly to 14 fs. 


VIII. RESULTS AND DISCUSSION 

The 7 -ray energies and mean lifetimes determined from 
the present study are listed in Table El and the ex¬ 
perimental 7 -ray spectra are shown in Figs. 11011201 with 
simulated line shapes superimposed. The fit range has 
been determined on a case-by-case basis, with the aim 
of maximizing the coverage while still allowing a sim¬ 
ple parametrization of the background as a constant or 
linear function plus a smoothed step function. In some 
cases, only upper limits on the lifetime could be ob¬ 
tained. The error bars represent our best estimate of 
the overall experimental uncertainty, including both the 
statistical uncertainty from the fit and systematic uncer¬ 
tainties. In most cases, the statistical uncertainty dom¬ 
inates. The systematic uncertainties were estimated in 
an ad-hoc fashion by varying the width of the a-particle 
gate, adopting different background parameterizations, 
widening/narrowing the fit range, and varying selected 
input parameters for the Monte Carlo simulation, such 
as the 7 -a angular-correlation function. The quality of 
the experimental data did not warrant a comprehensive 
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FIG. 10. Line-shape analysis of the a-gated, doppler-shifted, 
7 -ray line of the 2052 keV —> 451 keV transition in 23 Mg. The 
curves show simulated line shapes, which have been fitted to 
the experimental data (see text for details). The numbers in 
parentheses give the y 2 /dof of each fit. 



FIG. 11. The 2771 keV —> g.s. transition. 


and detailed sensitivity-analysis similar to the one re¬ 
ported in Ref. [22j]. We note that the non-observation 
of the 7787 keV —> 2051 keV transition in our experiment 
is consistent with its smaller branching ratio relative to 
the 7787 keV —> 451 keV transition and the higher 7 -ray 
background at lower energy. 

While we have succeeded in determining several, hith¬ 
erto unknown, lifetimes in 23 Mg, and placed relevant up¬ 
per limits on even more, our result for the 7787 keV state, 
r < 12 fs at the 95% C.L., is not sufficiently precise to 
settle the controversy surrounding the 22 Na(p, 7 ) rate. 
We note that our result is consistent with the mean life¬ 
time of r = 10 ± 3 fs reported in Ref. [llj], but also 
the mean lifetime of r = 2.4 +q ' 2 fs implied by the new 
resonance-strength determination of the Seattle group, 
wy = 5.7±o j meV @,i|, assuming B p = 0.037T0.07 [13], 

Our sensitivity has been limited mainly by the low 



FIG. 12. The 2908 keV —>- g.s. transition. 



FIG. 13. The 3798 keV -7 451 keV transition. 



FIG. 14. The 4353 keV —> g.s. transition. 
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TABLE II. Mean lifetimes, t, and 7 -ray energies, determined from line-shape fitting. States are labeled by their excitation 
energy in keV, as given in the most recent evaluation [3l|. Upper limits are given at the 95% C.L. 


Initial state 

Final state(s) 

E 1 (keV) 


T (fs) 

Present 

Literature 

2052 

451 

1601.4 ± 1.3 

104 ± 18 

80 ± 20 [321 

2771 

g.s. 

2769.7 ± 1.2 

98 ± 15 

155 ± 30 [321 

2908 

g.s. 

2906.4 ± 1.1 

15 ±3 

< 25 [32] 

3798 

451 

3343.9 ± 1.3 

41 ±6 

< 20 x 10 6 [31] 

4353 

g.s. 

4356 ± 2 

< 11 

< 20 x 10 6 [31] 

5287 

451 

4837 ± 3 

< 14 

5 ± 2 [121 

5453 

2052 

3401 ± 2 

< 15 


6236 

2052 

4186 ± 3 

< 40 


6375 

2052 

4322 ± 6 

< 45 


6899 

g.s. 

6907 ± 3 

< 10 


7444 

g.s. 

7443 ± 3 

< 14 


7493 

2052 

5442 ± 2 

< 20 


7787 

451 

7335 ± 2 

< 12 

10 ±3 [11] 



FIG. 15. The 5287 keV -7- 451 keV transition. 



Ey (keV) 


FIG. 16. The 5453 keV —> 2052 keV transition. 



FIG. 17. The 6375 keV —> 2052 keV transition. 



FIG. 18. The 6899 keV —> g.s. transition. 
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FIG. 19. The 7493 keV —> 2052 keV transition. 



FIG. 20. The 7787 keV -»• 451 keV transition. 


statistics collected on the 7787 keV —> 451 keV transition, 
which may be ascribed to a rather small cross section. 
Previous to our experiment the cross section for pop¬ 
ulating the 7787 keV state via ( 3 He,a) was unknown. 
Based on our data we estimate the cross section to be 
dcr/dflc.m. ~ 3-4 /jb/sr in the angular range covered in 
the present experiment (0 cm . > 159°) assuming that the 
7 rays are emitted isotropically. Our sensitivity has been 
further limited by the high-energy 7 -ray background re¬ 
sulting from reactions between the beam and the residual 
gas and water vapor that condensed on the surface of the 
cooled target, cf. Section HTT1 

It is also worth considering if the geometry of the de¬ 
tector setup can be improved. Fig. [21] (a) shows the sim¬ 
ulated line shape of the 7787 keV —> 451 keV transition 
with the current experimental setup assuming mean life¬ 
times of 0, 1, 5, and 15 fs. In Fig. [2T] (b) the same four life¬ 
times have been plotted, but with a modified setup: the 
angular acceptance of the A E-E telescope has been dou¬ 
bled while the angular acceptance of the HPGe detector 
has been halved (thus keeping the product f 2 Q fl 7 nearly 



FIG. 21. (a) Simulated line shape of the 7787 keV —> 451 keV 
transition with the current experimental setup assuming mean 
lifetimes of 0, 1, 5, and 15 fs. (b) Same, but with the accep¬ 
tance of the A E-E telescope doubled and the acceptance of 
the HPGe detector halved. 


constant). Based on the simulated line shapes it would 
appear that the lifetime sensitivity can be improved by 
simply reducing the acceptance of the HPGe detector, 
while increasing the acceptance of the A E-E telescope 
to maintain the same level of statistics. It must be re¬ 
membered, however, that the increased acceptance of the 
A E-E telescope will result in a worse excitation-energy 
resolution. The gates in the a-particle spectrum will have 
to be made wider and consequently more background will 
appear in the cc-gated 7 -ray spectra. Also, the increased 
acceptance of the A E-E telescope will increase the sensi¬ 
tivity to the angular distribution of the a-particles, which 
is usually not known. Thus, the preferred solution would 
be to keep the acceptance of the A E-E telescope fixed 
and employ several HPGe detectors, rather than just one, 
so that the HPGe detectors can be placed farther away, 
or employ a position sensitive HPGe detector such as the 
one described in Ref. 33|. This would, however, compli¬ 
cate both the experiment and the data analysis. 
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IX. CONCLUSION 

The 22 Na(p, 7 ) reaction plays an important role in clas¬ 
sical novae by limiting the production of the long-lived 7 - 
ray emitter 22 Na [lj. At the peak temperatures occurring 
in classical novae, the reaction rate is dominated by a sin¬ 
gle resonance due to an excited state in the compound nu¬ 
cleus, 23 Mg, at an excitation energy of 7787 keV. The two 
direct measurements of the resonance strength reported 
in the literature differ by a factor of 3 @|- The res¬ 

onance strength deduced indirectly from the properties 
of the 7787 keV state is consistent with the lower of the 
two reported direct measurements, but hinges on a single 
lifetime measurement with a 30% error bar El- Here we 
have reported on a new lifetime measurement performed 
at TRIUMF using the DSAM technique [l7l[l8|. We have 
developed a Monte Carlo simulation program to model 
the measured 7 -ray line shapes. Our program improves 
on a similar program previously used at TRIUMF ( 22 I ] 
by including the effects of multiple scattering. We have 
successfully determined several lifetimes in 23 Mg for the 
first time. Our result for the 7787 keV state, r < 12 fs 
at the 95% C.L., is unfortunately not sufficiently precise 
to settle the controversy surrounding the 22 Na(p, 7 ) rate. 
When combined with the proton-decay branching-ratio 
determination of Ref. [l3[, B p = 0.037 ± 0.007, and as¬ 
suming a spin of J = 7/2, our upper limit on the lifetime 
of the 7787 keV state gives a lower limit on the resonance 
strength of wy > 1.1 meV at the 95% C.L. 


ACKNOWLEDGMENTS 

We thank W. N. Lennard for helpful discussions on the 
properties of the 3 He-implanted Au target, and we thank 
E. Rand for supplying us with the GEANT4 simulation 
of the HPGe detector. OSK acknowledges support from 
the Villum Foundation. TRIUMF receives federal fund¬ 
ing via a contribution agreement through the National 
Research Council of Canada. 


13 


[1] J. Jose, M. Hernanz, and C. Iliadis, 

Nuclear Physics A 777 , 550 (2006), special Issue on 
Nuclear Astrophysics. 

[2] D. D. Clayton and F. Hoyle, Astrophysical Journal 187 , 

L101 (1974). 

[3] A. F. Iyudin, K. Bennett, H. Bloemen, R. Diehl, 

W. Hermsen, G. G. Lichti, D. Morris, J. Ryan, V. Schoen- 
felder, H. Steinle, A. Strong, M. Varendorff, and C. Win¬ 
kler, Astronomy and Astrophysics 300 , 422 (1995). 

[4] M. Hernanz, Bulletin of the Astronomical Society of India 
40 , 377 (2012). 

[5] D. C. Black, Geochimica et Cosmochimica Acta 36, 347 (1972) 

[6] S. Seuthe, C. Rolfs, U. Schroder, W. H. Schulte, E. So- 
morjai, H. P. Trautvetter, and F. B. Waanders, Nuclear 
Physics A 514 , 471 (1990). 

[7] F. Stegmiiller, C. Rolfs, S. Schmidt, W. H. Schulte, H. P. 
Trautvetter, and R. W. Kavanagh, Nuclear Physics A 
601 , 168 (1996). 

[8] A. L. Sallaska, C. Wrede, A. Garcia, D. W. Storm, 

T. A. D. Brown, C. Ruiz, K. A. Snover, D. F. Ottewell, 

L. Buchmann, C. Vockenhuber, D. A. Hutcheon, and 
J. A. Caggiano, Phys. Rev. Lett. 105 , 152501 (2010) 

[9] A. L. Sallaska, C. Wrede, A. Garcia, D. W. Storm, 

T. A. D. Brown, C. Ruiz, K. A. Snover, D. F. Ottewell, 

L. Buchmann, C. Vockenhuber, D. A. Hutcheon, J. A. 
Caggiano, and J. Jose, Phys. Rev. C 83, 034611 (2011) 

[10] A. Parikh, J. Jose, and G. Sala, 

AIP Advances 4 , 041002 (2014) 

[11] D. G. Jenkins, C. J. Lister, R. V. F. Janssens, T. L. 

Khoo, E. F. Moore, K. E. Rehm, B. Truett, A. H. 
Wuosmaa, M. Freer, B. R. Fulton, and J. Jose, 
Physical Review Letters 92, 031101 (2004) 

[12] D. G. Jenkins, M. Bouhelal, S. Courtin, M. Freer, B. R. 
Fulton, F. Haas, R. V. F. Janssens, T. L. Khoo, C. J. 

Lister, E. F. Moore, W. A. Richter, B. Truett, and A. H. 
Wuosmaa, Phys. Rev. C 87 , 064301 (2013) 

[13] A. Saastamoinen, L. Trache, A. Banu, M. A. Bent¬ 
ley, T. Davinson, J. C. Hardy, V. E. Iacob, M. Mc- 
Cleskey, B. T. Roeder, E. Simmons, G. Tabacaru, 

R. E. Tribble, P. J. Woods, and J. Aysto, 

Phys. Rev. C 83, 045808 (2011) 

[14] V. Tripathi, S. L. Tabor, A. Volya, S. N. Liddick, P. C. 
Bender, N. Larson, C. Prokop, S. Suchyta, P.-L. Tai, and 
J. M. VonMoss, Phys. Rev. Lett. Ill, 262501 (2013) 

[15] V. E. Iacob, Y. Zhai, T. Al-Abdullah, C. Fu, J. C. Hardy, 

N. Nica, H. 1. Park, G. Tabacaru, L. Trache, and R. E. 
Tribble, Phys. Rev. C 74 , 045810 (2006) 

[16] Y. Zhai, The Structure of Mg and Astrophysical Con¬ 
sequences, PhD thesis, Texas A&M University (2007). 

[17] D. Branford and I. Wright, 

Nuclear Instruments and Methods 106 , 437 (1973) 

[18] T. K. Alexander and J. S. Forster, Advances in Nuclear 
Physics 10 , 197 (1978). 

[19] J. Forster, D. Ward, H. Andrews, G. Ball, 

G. Costa, W. Davies, and I. Mitchell, 


Nuclear Instruments and Methods 136 , 349 (1976) 

[20] R. Kanungo, T. K. Alexander, A. N. Andreyev, G. C. 

Ball, R. S. Chakrawarthy, M. Chicoine, R. Church¬ 
man, B. Davids, J. S. Forster, S. Gujrathi, G. Hack- 
man, D. Howell, J. R. Leslie, A. C. Morton, S. Mythili, 

C. J. Pearson, J. J. Ressler, C. Ruiz, H. Savajols, 

M. A. Schumaker, I. Tanihata, P. Walden, and S. Yen, 

Phys. Rev. C 74 , 045803 (2006) 

[21] S. Mythili, B. Davids, T. K. Alexander, G. C. Ball, 

M. Chicoine, R. S. Chakrawarthy, R. Churchman, 

J. S. Forster, S. Gujrathi, G. Hackman, D. Howell, 

R. Kanungo, J. R. Leslie, E. Padilla, C. J. Pear¬ 
son, C. Ruiz, G. Ruprecht, M. A. Schumaker, I. Tani¬ 
hata, C. Vockenhuber, P. Walden, and S. Yen, 

Phys. Rev. C 77 , 035803 (2008) 

[22] N. Galinski, S. K. L. Sjue, G. C. Ball, D. S. Cross, 

B. Davids, H. A1 Falou, A. B. Garnsworthy, G. Hack- 
man, U. Hager, D. A. Howell, M. Jones, R. Kanungo, 

R. Kshetri, K. G. Leach, J. R. Leslie, M. Moukaddam, 

J. N. Orce, E. T. Rand, C. Ruiz, G. Ruprecht, M. A. 

Schumaker, C. E. Svensson, S. Triambak, and C. D. 

Unsworth, Phys. Rev. C 90 , 035803 (2014) 

[23] J. F. Ziegler, J. P. Biersack, and M. D. Ziegler, SRIM - 
The Stopping and Range of Ions in Matter, 5th ed. (SRIM 
Co., United States of America, 2008). 

[24] M. S. Basunia, Nuclear Data Sheets 114 , 1189 (2013) 

[25] R. Anne, J. Herault, R. Bimbot, H. Gau- 

vin, G. Bastin, and F. Hubert, 

Nuclear Instruments and Methods in Physics Research Section B: B< 

[26] D. Guillemaud-Mueller, M. O. Lam- 

pert, D. Pons, and M. Langevin, 

Nuclear Science, IEEE Transactions on 33, 343 (1986) 

[27] R. Brun and F. Rademakers, 

Nuclear Instruments and Methods in Physics Research Section A: A< 
new Computing Techniques in Physics Research V. 

[28] C. Svensson and A. Garnsworthy, 

Hyperfine Interactions 225, 127 (2014) 

[29] T. Alexander, G. Ball, W. Davies, and I. Mitchell, 

Journal of Nuclear Materials 96, 51 (1981) 

[30] S. Baker and R. D. Cousins, 

Nuclear Instruments and Methods in Physics Research 221 , 437 (19< 

[31] R. B. Firestone, Nuclear Data Sheets 108 , 1 (2007) 

[32] R. Engmann, F. Brandolini, and I. Mauritzson, 

Nuclear Physics A 171 , 418 (1971) 

[33] M. Schumaker, G. Hackman, C. Pearson, C. Svens¬ 
son, C. Andreoiu, A. Andreyev, R. Austin, G. Ball, 

D. Bandyopadhyay, A. Boston, R. Chakrawarthy, 

R. Churchman, T. Drake, P. Finlay, P. Garrett, 

G. Grinyer, B. Hyland, B. Jones, R. Maharaj, A. Mor¬ 
ton, A. Phillips, F. Sarazin, H. Scraggs, M. Smith, 

J. Valiente-Dobon, J. Waddington, and L. Watters, 

Nuclear Instruments and Methods in Physics Research Section A: A< 





